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ABSTRACT 

Observations of the most distant bright quasars imply that billion solar mass supermassive black holes 
(SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy for- 
mation scenario such fast growth implies large gas densities providing sustained accretion at critical or super- 
critical rates onto an initial black hole seed. It has been a long standing question whether and how such high 
black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new 
MassiveBlack cosmological hydrodynamic simulation covering a volume (0.75 Gpc) 3 appropriate for studying 
the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies 
produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensu- 
rate with the existence of ~ 10 9 M black holes as early as z ~ 7. We find that under these conditions quasar 
feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench 
the accretion until the host galaxy reaches Mh a i ~ 10 12 M . This cold-flow driven scenario for the formation 
of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy 
mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs. 
Subject headings: quasars: general — galaxies: formation — galaxies: active — galaxies: evolution — cos- 
mology: theory — hydrodynamics 



1. INTRODUCTION 

It is now well established that the properties of supermas- 
sive black holes (SMBH) found at the centers of galaxies 
today are tightly coupled to those of their hosts implying a 
strong link between black hole and galaxy formation. The 
strongest direct constraint on the high-redshift evolution of 
SMBHs comes from the observations of the luminous quasars 
at z ~ 6 in the Sl oan Digital Sky Survey (SDSS) (Fa net al.l 
2006; Jiang et al. 2009) and even more recently at z = 7 
(Mortloc k et al .11201 1|) . Although rare (the comoving space 
density of z ~ 6 quasars is roughly n ~ a few Gpc -3 ) the in- 
ferred hole masses of these quasars are in excess of 10 9 M 
comparable to the masses of the most massive black holes in 
the Universe today. The origin of these massive black hole 
seed and the physical conditions that allow early growth to 
supermassive black holes remain a challenging problem. 

In order to have had sufficient time to build up via gas accre- 
tion and BH mergers (resulting from the hierarchical merging 
of their host halos) the first 'seed' black holes must have ap- 
peared at early epoch, z > 10. The origin and nature of this 
seed population remain uncertain. Two distinct populations 
of seed masses, in the range of 100 — 10 5 M have been pro- 
posed: the small mass seeds are usually thought to be the rem- 
nants of the first gener ation of PopIII stars formed of metal- 
free gas at z ~ 20 - 30 (Bromm etaLllT99 9: Abel et al. 2000; 
Naka mura & Umem ura 2001; Yoshidaetal. 2003; Ga o et al.l 
2006, e.g.), while the large seeds form in direct dy nami- 
cal collapse in metal- free galaxies (Koushiap pas et al.l 12004: 
iBegelman et al.l I2006L although see also Mayer et al. 2010 
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for direct collapse into a massive blackhole in metal enriched 
regime). 

Growing the seeds to 1O 9 M in less than a billion years 
requires extremely large accretion rates - as mergers between 
black holes are too rare and too inefficient for significant 
growth. For a black hole accreting at the critical Edding- 
ton accretion rate its luminosity L Edd = {AirGcm p ) / a T M BH = 

r]MEddC 2 (where G, c, m p and a? are the gravitational con- 
stant, speed of light, proton mass and Thomson cross sec- 
tion, and r] is the standard accretion efficiency) implies an ex- 
ponential growth at the characteristic Eddington timescales, 
fedd = 450^/(1 - rj) Myr, such that M BH = M SQQd e t/tEdd . For a 
seed mass ranging from M see d ~ 100 — 1O 5 M this requires 
10—17 e-foldings to reach Mbh ~ 10 9 M . The crucial ques- 
tion (for any given seed model) is then where (if at all, in 
which kind of halos) and how (at what gas inflow rates) such 
vigorous accretion can be sustained at these early times. 

As bright quasars are likely to occur in extremely rare high- 
density peaks in the early universe, large computational vol- 
umes are needed to study them. Here we use a new large 
cosmological Smooth Particle hydrodynamics (SPH) simula- 
tion, Massive Black (covering a volume of [0.75 Gpc] 3 ) with 
sufficiently high resolution (over 65 billion particles) to be 
able to include tested prescriptions for star formation, black 
hole accretion and associated feedback processes to investi- 
gate whether and if such objects may be formed within our 
standard structure formation models. Crucially, our Massive - 
Black simulation is of sufficiently high-resolution to allow us 
to follow the mass distribution in the inner regions of galax- 
ies and hence model star formation and black hole growth di- 
rectly and self-consistently while still evolving a close to Gi- 
gaparsec scale region. It therefore provides a unique frame- 
work to study the formation of the first quasars. 



2. METHODOLOGY 
2.1. Simulation run: "Massive Black' 
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Fig. 1 . — The cosmological mass distribution in our the simulation volume at z = 5. The projected gas density over the whole volume ('unwrapped' into 2D) 
is shown in the large scale image. The two overlaid panels show successive zoom-ins by factor of 10, center on the region where the most massive black hole is 
found. 



Our new simulation has been performed with the cos- 
mological TreePM- Smooth Particle Hydrodynamics (SPH) 
code P -GADGET, a hy brid version of the parallel code GAD- 
GET2 (Springel 2005) which has been extensively modified 
and upgraded to run on the new generation of Petaflop scale 
supercomputers (e.g. machines like the upcoming Blue Waters 
at NCSA). The major improvement over previous versions of 
GADGET is in the use of threads in both the gravity and SPH 
part of the code which allows the effective use of multi core 
(currently 8-12 cores per node) processors combined with an 
optimum number of MPI task per node. Here we present ini- 
tial results from the largest cosmological smooth particle hy- 
drodynamic simulation to date which was run on 10 5 cores 
corresponding to the entire Cray-XT5 "Kraken" at NICS. The 
MassiveBlack simulation contains N part = 2 x 3200 3 = 65.5 
billion particles in a volume of 533 Mpc//z on a side with 
a gravitational smoothing length e = 5.5kpc//z in comoving 
units). The gas and dark matter particle masses are m g = 
5.7 x 1O 7 M and m DM = 2.8 x 1O 8 M respectively. This 
run contains gravity and hydrodynamics but als o extra physics 
(subgrid modeling) for star formation (Spring el & Hernquistl 
2003), black holes and associated feedback processes. The 
simulation has currently been run from z=159toz = 4.75 
(beyond our original target redshift of z = 6). For this mas- 
sive calculation it is currently prohibitive to push it to z = as 
this would require an unreasonable amount of computational 
time on the world's current fastest supercomputers. The sim- 
ulated redshift range probes early structure formation and the 
emergence of the first galaxies and quasars. 

2.2. Black hole Accretion and Feedback Model 

The prescription for accretion and associated feed- 
back from ma ss ive b la ck holes has been developed by 
iDi Matteo et~aTI (120051) : ISpringel et all (120051) . Detailed 
studies of this implementation in cosmological simula- 
tions and associated predictions ( Sij acki et al.ll2007l ; iLi et al] 



2007; Di Matteo et al. 2008; Croft etal. 2009; Siiacki et al. 
2009; Colbert & di Matteo 2008]; iDegraf et al] 1201 a 12011 bl: 
Booth & Schaye 2011) have shown that it can reproduce 
all the basic properties of black hole growth, the observed 
Mbh — cr, relation (Di Matteo et al. 2008), the quasar lumi- 
nosity function (Degraf et al. 2010) and its evolution as well 
as the spatial clustering of quasars dDeGraf et al] l2011cl) . 
In a nutshell our blac k hole accretion and feedback model 
(IDi Matteo et al .120081) consists of representing black holes by 
collisionless particles that grow in mass (from an initial seed 
black hole) by accreting gas in their environments. A frac- 
tion of the radiative energy released by the accreted material 
is assumed to couple thermally to nearby gas and influence 
its motion and thermodynamic state (typically referred to as 
BH feedback). Our underlying assumption is that the large- 
scale feeding of galactic nuclei with gas (which is resolved 
in our simulations) is ultimately the critical process that de- 
termines the growt h of massive black holes and the peak of 
the q uasar phase (IDi Matteo et al. 1 120081: IDegraf et al.ll201QL 
201 lb). The model, therefore needs to be viewed in the con- 
text of cosmological growth of black holes and not detailed 
accretion physics. While a more detailed treatment of this 
is certainly de sirable and begins to be possible for i ndividual 
galaxies dKim et al.ll2QTTt iHopkins & Ouataertll201Q[ e.g.) it 
is still infeasible for cosmological simulations that seek to fol- 
low whole populations of galaxies and their BHs. 

We introduce collisionless 'sink' particles in the simula- 
tions to model black holes at the centers of forming minihalos. 
In order to achieve this we keep track of the formation of mini- 
halos by running a friends-of-friends (FOF) group finder on 
the fly. The group finder is run on sufficiently closely spaced 
intervals to identify the newly collapsing halos in which we 
place a black hole seed of fixed mass, M— 10 5 /z -1 M (if they 
do not already contain a BH). In practice most of the halos of 
this mass and their black holes are formed between z = 15 
and z = 30. The black hole particle then grows in mass via 
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Fig. 2. — Snapshots of the evolution of the cold-flow-fed massive black holes. The images visualise the projected gas distribution color coded by temperature 
around three example quasars (one in each row) across five different redshifts (labelled on each of the five panels; left to right). The projected density ranges from 
~ 10 -2 to ~ 10 2 hM© pc -2 and the temperature from ~ 10 4 (blue colors) to 10 8 K (red colors). The quasar positions are indicated by the green circles and the 
virial radius of each halo by the blue circles. The images show the typical structure of cold streams that penetrate the halo all the way into the central regions of 
galaxies. At first (z ^ 7.5) the gas is cold and the black hole is still of relatively low mass, but then below this redshift the black hole growth exponentiates as 
increasing amounts of cold gas is fed into the central regions. Black hole feedback heats the gas, but does not distrupt the cold streams until z ^ 6. 



accretion of surrounding gas according to M B h = / 2+v 2W2 
(where p and c s are the density and sound speed of the hot 
and cold phase of the ISM g as which whe n taken into account 
appropriately as in Pelupess y et al.l (|2007|) - this eliminates the 
need for a correction factor a previously introduced - and vbh 
is the velocity of the black hole relative to the gas) and by 
merging with other black holes. We limit the accretion rate to 
Eddington (or a few tim es Eddington). A s imilar Bondi model 
has also been used by I Johnson & Bromrnl (120071) to study the 
growth of the first massive PopII black holes remnants. 

The radiated luminosity, Lr, from the black hole is related to 
the accretion rate, Mrh as Lr = e r (Mbh x c 2 ) , where we take 
the standard mean value e r = 0.1. Some coupling between 
the liberated luminosity and the surrounding gas is expected: 
in the simulation 5% of the luminosity is (isotropically) de- 
posited as thermal energy in the local black hole kernel, pro - 
viding some form of feedback energy (Di Mat teo et"aT1l2005l) . 
This model of AGN feedback as isotropic thermal coupling 
to the surrounding gas, albeit simple, is a reasonable approx- 
imation to any physical mechanism which leads to a shock 
front which isotropizes and becomes well mixed over physi- 
cal scales smaller than those relevant in our simulations and 
on timescales smaller than the dynamical timescales of the 
halos (Di Ma tteo et al] [2QQ8: Hopkins et al. 200^). Two black 
hole particles merge if they come within the spatial resolu- 
tion (i.e. within the local SPH smoothing length) with relative 



speed below the local sound speed. 

We note that at leas t two ind ep enden t 

groups dBooth & Schavel 120091: iJohansson et~aTl 12008b 
now have also adopted the same modeling for black hole 
accretion, feedback and BH mergers in the context of hy- 
drodynamic simulations. These independe nt works, and in 
partic ular, the cosmological simulations by (JBooth & Schavel 
|2009|) (part of the OWL program) have allowed to indepen- 
dently explore th e parameter space of the reference model of 
(Di Matte o et al.l [2008), as well as variations of our model 
prescriptions. This large body of already existing work and 
investigations make this particular model a good choice for 
more detailed studies of the growth of the first quasars which 
is the subject we focus on here. In associated publications we 
show that the MassiveBlack quasars are fully consistent with 
all the fundamental statistical constraints for the observed 
populations of high redshift quasars, and in particular with 
the observed luminosity functio ns of q uasars and the high 
redshift clustering (DeGraf et al.ll201 la|) and basic properties 
of quasar hosts (|Khandai et al.ll201 1|) . To produce the black 
hole history from the simulation data we rely on SQL 
databases developed by lLopez et al.l (l201l|) . 

3. RESULTS 

The cosmological gas density distribution in the full volume 
of the MassiveBlack is shown in the large scale image of Fig- 
ure 1 . The large panel shows the whole of the 3D simulation 
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Fig. 3. — The black hole mass and the black hole accretion rate versus redshift are shown in the top and bottom panel on the left. The corresponding 
densities and temperatures of the accreting gas are shown in the top and bottom panels on the right. The lines represent the three black holes shown in Fig. 2 
the shaded orange band indicates the full range of properties encompassed by our most massive black hole sample. 
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volume "unwrapped" into a 2D image slice (|Feng et al .11201 1|) 
at z = 5. At these large scales the density distribution of the 
universe appears fairly uniform. The resolution however is 
sufficiently fine to make it possible to zoom into increasingly 
smaller regions and search for massive black holes that have 
experienced significant growth. Superimposed onto the full 
scale image we show a zoomed region (scale of lOOMpc h~ l 
and lOMpc h~ l on a side from left to right, respectively) 
around one of the largest black holes/quasars at this time. 
On these scales the images show the typical filaments that 
compose the cosmic web and in particular how the first mas- 
sive quasars (the most massive black hole, at the center has 
a mass M B h ~ 3 x IO 9 M ) form at the same type of inter- 
section/nodes of filaments that are the expected locations of 
rare (massive) dark matter halos. Remarkably we do find ten 
black holes in the volume that have grown to about a billion 
solar masses by z ~ 6 or earlier (and many other of smaller 
masses). 

Figure 2 shows the environment and its evolution (five time- 
lapses from about z~9toz~5)of three examples from 
within this sample. Their detailed mass assembly history is 
shown in Fig. 3. The panels in Figure 2 show the evolu- 
tion of the gas density color coded by temperature. These 
objects are found to be continuously fed by intense streams 
of high density gas (consistent with the cold- accretion pic- 
ture f or the growth of galaxies at intermediate/high redshift 
by (iDekel & Birnboimll2006l e.g.). During these times Ed- 
dington accretion is attained and sustained. By showing the 
temperature of the gas, the images clearly make visible an ex- 
panding 'bubble' (emerging from about z ^ 6.5 — 7) of hot 
gas (red colors) around the central quasars (whose positions 
are indicated by green circles). This bubble created by the BH 
feedback is more or less confined within the halo (the virial 
radius of the halo is shown by the blue circles in Fig. 2) for 
Z ;$ 6. Below this redshift, and rather abruptly, the energy re- 
leased by the quasars heats and expels the gas as a wind well 
beyond the halo. Black hole growth has now become self- 
regulated (see also Figure 3, Eddington rates are reached only 
sporadically). Although the effects of quasar feedback in our 
model have been studied in detail previously (Di Matteo et al. 



2005, 2008) what is remarkable here is that even though feed- 
back energy consistently heats the gas within and eventually 
beyond the scale of halos, it does not do much to the streams 
of in-flowing cold gas. There the gas density is so high (e.g. 
Fig. 3) that the gas cannot be stopped because, it is too dif- 
ficult to couple enough feedback energy to it to disrupt the 
flow. The streams get somewhat (albeit not completely) dis- 
rupted only at z ^ 6. Before this happens billion solar mass 
black holes are already assembled (see Fig. 3), a process that 
takes a few hundred million years. 

Sustained phases of Eddington accretion onto these black 
holes start as early as z ~ 9 — 10 and go on uninterrupted 
until z ~ 6 — 7, leading to BH masses of the order of 
IO 9 M in the first massive halos of roughly Mh a i ~ IO 12 M 
at z ~ 6 — 7 (rare 3 — 4a peaks of the density distribu- 
tion). We have indeed a few objects (the range of the 10 
most massive is shown by the orange areas in Fig. 3) that 
reach these high ma sses already at z ~ 7 consistent with 
(Mortloc k et al.ll20TT|) . Interestingly, this process of BH as- 
sembly is apparently facilited by the "cold flows" picture of 
galaxy formation which has revolutionized our understand- 
ing of galaxy assembly below the threshold dark matter mass 
of Af h aio ~ IO 12 M^ dDekel & BirnboTml 120061: IDekel et al.l 
120091: JKeres et al.ll2005l 120091) . With MassiveBlack we are 
able to trace the formation of the first, rare massive halos, 
those which are mostly assemble d by high density (high red- 
shift) cold streams (IDekel et al.ll2009h . We find that these 
same streams easily penetrate all the way into central re- 
gions of galaxies even in the presence of strong feedback. We 
find however that once Mh a i ~ IO 12 M and halos enter the 
regime where they are shock-heated (Dekel et al. 2009) the 
BH growth becomes finally self-regulated (see accretion rate 
evolution in Fig. 3). The temperature of the accreting gas also 
is raised well above the virial temperature, T yiv ~ IO 7 K, ren- 
dering some of the gas unbound; Fig. 3). This the point the 
black hole masses level off at a few 1O 9 M (Fig. 3). Even 
though at z < 6 in massive halos the outflows clearly affect the 
incoming gas it is not clear they fully d istrupt the cold flows. 
According to IDekel & Birnboiml (I2006L e.g), albeit based on 
somewhat lower redshift than what we probe, cold streams 
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Fig. 4. — The temperature of the gas that is eventually accreted onto the 
black hole is shown as a function of its position within the host halo (r is 
plotted in physical units). The plot shows a typical example for one of the 
host halos in Fig. 3. The points with different symbols representing the 
groups of 32 nearest neighbours at different redshifts (in the colors indicated 
by the legend) all traced back to z = 8.5. The gray hatched region show the 
gravitational softening lengths 



should still not be fully distrupted in this rare massive objects 
even if a dilute hot medium forms. In future work we will 
address this issue in our simulations in some more details. 

We note that the black holes occupying these rare massive 
halos hardly undergo any mergers, indicative that their hosts 
are not undergoing major galaxies mergers. In the inset in Fig. 
3 we show the histories of all black holes that merge into the 
main progenitor, which is typical of what we see: a few and 
typically minor mergers occur. At the redshifts relevant for 
the first quasars, however, galaxy formation is very different: 
mergers are still extremely rare events and halos have not ye t 
assembled above their shock heating scale (De kel et al1l2 009). 
so that quasars are fueled directly by cold streams/flows. 

To illustrate directly the origin of the gas fueling the first 
quasars we track in the simulation the prior history of parti- 
cles that end up in the vicinity of the black hole, and hence 
contribute to the accretion. Figure 4 shows one example of 
the temperature as a function of radius of particles that have 
participated into the accretion onto the BH between z = 6 and 
z = 8.5 (within the smoothing length of the black hole parti- 
cle). The plot shows that particles always remain cold as they 
enter the virial radius well into the ten of kpc region in-fact all 
the way into the region of from which accretion onto the black 
hole occurs. Here some of the gas is heated, as expected, by 
BH feedback. The temperature of the gas feeding the black 



hole is fully commensurate with that of cold flows and always 
remain s below r gas ~ iq5-5-6 k j u s ^ as predicted in g alaxy for- 
mation |DekeI&iEnboi3 (l2QQ6h : iKeres et al.l (12005b . In the 
shock-heated regime the temperature of the gas would rise to 
r v i r as it enters the halo. However this is not seen. 



4. CONCLUSIONS 

With our new large cosmological simulation MassiveBlack 
we show that the short timescale associated with infall via 
cold flows and the short cooling timescales in cold radial 
streams that penetrate the halo render the flow into the cen- 
tral regions unstoppable by feedback allowing it to easily sus- 
tain BH growth at the Eddington rates to build up the required 
BH masses by z = 6 — 7. One consequence of this scenario 
is that BH masses at these redshifts are expected to show 
deviations from the local M B h — & relation. BH masses as- 
sembled faster (most growth occurs over few hundred mil- 
lion years) than the stellar spheroid (assembled over Gyrs 
timescales). Black hole masses larger than that inferred from 
the local M B h — cr relation have infact been suggested for the 
first quasars. Stream fed accretion will still be relevant over 
the peak of the quasar phase but mergers (as merger rates peak 
closer to those redshift) will become an increasingly major 
player in their growth and formation. We speculate however 
that most of the growth of a quasar's mass is likely to always 
occur before the shock heating scale of an halo is r eached 
much like most of its star formation rate (De kel et al.1 12009). 
We will investigate this further in our large volume in future 
work. Our scenario is somewh a t simil ar to that proposed by 
iMaver et al.l (1201 Oh or iLi et al.l (120071) yet, crucially it does 
not rely on a major merger to induce the strong gas inflows 
but points to a more common origin for them particularly at 
these redshifts (which we could find b y virtue of haying a 
large volume in our simulation, see also ISijacki etail (2009) 
who followed the build-up of a single hihigh-z quasar). As 
we have shown, relaxing the constraint for massive mergers 
makes it plausible to attain quite commonly la rge black hole 
masse s as high as z = 7 commensurate with Mo rtlock et al.l 
d20TTh . 
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